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We address the issue of internal modes of a kink of a discrete sine-Gordon equation. The main 
point of the present study is to elucidate how the antisymmetric internal mode frequency dependence 
enters the quasicontinuum spectrum of nonlocalized waves. We analyze the internal frequency 
dependencies as functions of both the number of cites and discreteness parameter and explain the 
origin of spectrum peculiarity which arises after the frequency dependence of antisymmetric mode 
returns back to the continuous spectrum at some nonzero value of the intersite coupling. 
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The existence of internal modes is a well known pe- 
culiarity of nonintegrable systems 0] : If one linearizes a 
nonintegrable equation above the static nonlinear solu- 
tion, then the effective potential of the linear equation 
obtained appears to be not reflectionless in contrast to 
the integrable case. The particular example is the ex- 
istence of internal (shape) modes of kinks in the well 
known discrete sine-Gordon equation (DSGE). This set 
of differential-difference equations reads as (in normal- 
ized dimensionless units): 

u n + sinu n + A (u n - u n -i) - (u n +i - u n ) = 0, (1) 

where u n (t) is the field variable, which can have a multi- 
tude of physical meanings 0, A is a coupling parameter, 
dot means the time derivative and index n numerates the 
ID chain sites. The features and behaviour of kink inter- 
nal modes for DSGE, as well as for more general types of 
the so-called Frenkel-Kontorova (FK) model, have been 
already studied in very details 0. [J. |j. . The spectrum 
of linear waves of DSGE around a kink contains either 
two or one localized mode, depending on the value of 
parameter A. (We consider the stable kink centered be- 
tween the chain sites.) The frequencies of these modes 
lie inside the spectrum gap; the linear spectrum itself is 
given by 



1 + 4A sin 2 



2 ' 



(2) 



where uj is the frequency of linear waves and k is the wave 
number. The lowest (gap edge) frequency is u> = 1 (in 
the renormalized units). 

However, up to our knowledge, notwithstanding the 
great quantity of results none so far has concentrated on 
the way an internal modes ceases to exist: Only the fact 
of existence or nonexistence of that mode has been the 
subject of interest 0, S H S @. In this paper we focus 
on the effect entailed by the detachment of the internal 
mode frequency from the spectrum. So, the question 
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arises: How the internal modes, which split off from the 
spectrum at some nonzero value of coupling parameter 
(or some other effective parameter altering the system 
state), affect the spectrum itself and how they behave 
while still being inside the spectrum. 

To begin with, we note that the analytical approximate 
solution for a static kink of DSGE JTJ can be found either 
in the so-called anticontinuum limit (extremely small val- 
ues of A) as a series in powers of A 4] , or in the opposite 
case of strong coupling (large A) , when the discrete kink 
solution acquires the form of that in the continuous SGE 
with small corrections due to discreteness (see e.g. jy|). 
Let us now summarize the results on the existence of in- 
ternal modes in the DSGE containing a single kink. For 
small A only the lowest symmetric Pierls-Nabarro (PN) 
mode, associated with the kink oscillations in the PN 
relief, has the frequency inside the gap: it corresponds 
to the translational kink mode of continuous SGE, ac- 
tivated due to discreteness. This mode remains the in- 
ternal mode for the whole interval of kink stability. The 
roperties of this mode are well understood and studied 
|, and we shall be mainly concerned with the other 
internal mode. For larger (but still weak) A there exists 
the "critical" point A^, where one more mode dependence 
detaches from the continuous spectrum: this mode cor- 
responds to the oscillation of kink width. (More compli- 
cated on-site FK-type potentials may have a larger va- 
riety of localized internal eigenmodes Q or, for the case 
of more than one kink in the chain, the DSGE possesses 
a larger number of internal modes as well 0.) The ex- 
istence of the second internal mode for large A obviousl 
matches the criterium given by Kivshar et al. in Ref.Qj 
for nearly integrable SGE (see also Ref.[5| for details). 

The effects we shall discuss can be readily seen in Fig^ 
As the internal mode frequency goes back for small A to 
the quasicontinuous spectrum, it brings about a conspic- 
uous change in the behaviour of higher modes. However 
only the dependencies of odd modes (the mode number 
corresponds to the number of eigenfunction nodes) "feel" 
the return of the localized first antisymmetric mode to 
the spectrum. The modes with the even number of eigen- 
function nodes do not react to this return. So, we can in- 
fer that only the modes having the same symmetry as the 
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FIG. 1: Numerical results for DSGE containing 250 sites 
with a kink centered in the middle of the chain between the 
sites, a) The frequency dependencies for several lowest modes 
on the value of A (the dependence for the lowest PN-mode is 
not presented) ; the numbers correspond to the mode number, 
i.e. to the number of nodes of corresponding eigenfunction. 
The curve marked "1" is the dependence for localized an- 
tisymmetric mode. The inset shows the magnified region in 
the close vicinity of the detachment value Xd (the bandgap 
threshold marked by the dashed line) . The bottom panels 
show the behaviour of b) du/dX and c) d 2 w/dA 2 for three 
lowest antisymmetric modes (1,3 and 5) . 



mode, which detaches from the spectrum, experience pro- 
nounced change (which can be also seen in the behaviour 
of the derivatives, Fig|IjD,c). This change "propagates" 
inwards the spectrum becoming less and less abrupt with 
the increase of mode number. In addition, as can be seen 
in the inset of Fig^i, the dependence of the first anti- 
symmetric mode indeed enters the spectrum: it crosses 
the band edge frequency u) = 1 at the point A^ ~ 0.26 
(of course, it crosses only the band threshold but not 
any of the frequency dependencies). For smaller A the 
dependence of former antisymmetric localized mode be- 
longs to the spectrum of usual delocalized waves, tending 
to coalesce with the dependence of second (symmetric) 
mode. Also we note another interesting feature of eigen- 
frequency dependencies. For large A, i.e. when we are ap- 
proaching the continuum limit, the dependencies for odd 
modes with numbers 2n + 1 tend to coalesce with the de- 
pendencies for lower even modes, i.e. with those having 
numbers 2n. However the detachment of one mode re- 
sults in another tendency: after the splitting, for A close 
to zero, the dependencies for odd modes tend to coalesce 
with those of higher even modes, i.e. the frequencies of 
modes with numbers 2n + 1 come closer and closer to 
those for the modes with numbers 2n + 2. The inflec- 



tion of the odd modes dependencies and peculiarities of 
their derivatives in the vicinity of Xd merely indicate this 
"change" in the tendency for coalescence. 

To delineate and compare the behaviour of internal 
modes we first study some illustrative examples when we 
have only a few pendulums with linear bonds. 

Small number of pendulums with linear bonds. For the 
simplest case of two pendulums the dynamical equations 
read as: 



Mi, 2 + A(Mi i2 - "2,i) + sinui^ = . 



(3) 



Linearizing Eqs.© around the ground state solution, 
U\ = U2 = 0, we arrive at the eigenfrequency depen- 
dencies: uJq — 1 , lo\ — 1 + 2A. The higher frequency 
corresponds to the antiphase oscillations, the lower one - 
to the inphase motion. Now let us study the spectrum of 
linear oscillations above the inhomogeneous static solu- 
tion of Eq.© (analogue of a kink). First we have to find 
the static distribution for small values of A. We suppose 
that this distribution has a form: u\ = it , U2 = 2tt — u°, 
u° = aX + /3A 2 + 0(X 3 ), and we should determine the cor- 
responding constants, a and (3, in the dependence u°(A). 
Expanding the sine in Taylor series and equalizing coef- 
ficients at the same powers of A we find: 



i° = 2ttA - 4ttA 2 



(4) 



Then we linearize the initial set of dynamical equations 
© with respect to small quantities v^. u\ = u° + vie 1LOt , 
U2 = 2ir — u° + V2e luJt . The consistency condition gives 
the spectral dependencies: for the lower inphase mode 
luq = cos it , that is the analogue of PN mode in a large 
system, and for the antiphase mode uj 2 — cos u° + 2X — 
the analogue of the second antisymmetric internal mode. 
The analytical form for the frequencies, which is valid for 
the small A values, is as follows: 



uil = 1 - 2tt 2 A 2 



LJ? = UJn + 2A . 



(5) 



The frequency dependence wi(A) grows linearly with A 
for weak coupling parameter values and stays inside the 
band, but then, as A becomes larger, this growth slows 
down showing an ultimate tendency of this mode to get 
inside the gap below. The value of A for which this de- 
pendence crosses the band threshold can be roughly esti- 
mated as Xd ~ 1/7T 2 . One also notes that the frequency 
of PN-mode, loo, drops to zero for u° = ir/2. This man- 
ifests the well known fact that for the finite size chain 
there occurs an instability when the characteristic spa- 
tial scale of a kink overgrows the size of a system (see 
e.g. Ref.0, 8]). For a large system the PN-mode softens 
at some critical value A c ~ L 2 , where L 3> 1 is the sys- 
tem size 0. For two pendulums, setting u\ — it/ 2 — U, 
uPj = 3tt/2 + U, U <C 1, one finds the frequency depen- 
dencies near the critical point: w 2 (A c — A)/2A 2 , with 
A c = l/V, and ui\ = uiq + 2X. We see that at the criti- 
cal point u) 2 (X c ) — 2/tt < 1, i.e. the frequency remains 
inside the gap. Thus we expect that this internal mode 
dependence for a big system does not return back to the 



3 



spectrum for large values of A and continues to be the 
internal mode for the whole interval of kink stability ex- 
cept for the extremely weak coupling. The splitting of 
this mode dependence from the continuum spectrum is 
of the order of A -1 for large values of coupling parameter 
00. 

The systems containing more sites enable us to study 
the " 'degree of localization" since we now have the outer 
and inner (core) sites. Let us mark the sites in the kink 
core (the middle of a chain) by numbers "1" and "-1", 
their neighbor sites by "2" and "-2" etc. For the case of 
four pendulums the "kink" solution for small A is given 
by: 

u° 2 « 2ttA 2 + . . . , u\ « 2ttA - 8ttA 2 + . . . , (6) 

and the obvious solution symmetry, u-i — 2tt — Ui, gives 
the field distribution for the remaining sites. In this limit 
the spectral dependencies for two lowest modes are: 

ujI rs 1 - tt 2 A 2 , (7a) 
w^l + (2-^)A-y(2-V2)A 2 . (7b) 

Putting uji{\) = 1 in Eq.j7E|) one gets the rough estimate 
for the point where this frequency dependence enters the 
" quasicontinuum" spectrum: Xd w 2/ir 2 . In the same 
way one could deal with the systems of bigger size (see 
below). 

The spectrums and amplitude ratio dependencies on 
A for four and six pendulums are given in Fig[21 As the 
number of sites grows, the spectrum peculiarities, seen in 
Fig'H become more and more pronounced. In the panels 
(b) and (d) of Fig[3]we present the dependencies for the 
ratios of the oscillation amplitudes for different pairs of 
neighbor sites, Vi+i/vi, for PN (zeroth) mode and first 
antisymmetric mode. Both modes display the tendency 
for localization in the kink core with the growth of A, but 
coming closer to the critical point A c they become local- 
ized worse because of the large kink width. Compared 
with the first antisymmetric mode, the PN-mode appears 
to be localized better. These dependencies have a min- 
ima for some intermediate value of coupling, where the 
inhomogeneity brought by the kink is pronounced well 
enough to sustain the modes localization and, from the 
other side, the kink wight is still not large. 

Large systems. Now suppose that our system has 27V 
sites. The static distribution for the kink is again given 
by = 2ttA + C(A 2 ), u° = 0(X n ). Then we substi- 
tute u n {t) — u° + v n e luJt , and linearize the system ob- 
tained with respect to v n noting that for the antisym- 
metric modes the symmetry is v n = —v- n . After that 
one gains the following linear system: 

(u 2 — 3A — cos u\)v\ + \v 2 = , 
(w 2 -2\-CD8U°)v n + \(v n _ 1 +v n+1 ) = 0, (8) 
(u> 2 — A — cos u° n )vn + \vn-i = . 




0,0 0,2 y 0,4 0,6 



CO 

1,5 

C)i.o 





""3 








1 


6 pendulums 





0,0 0,25 ^ 0,75 




0,0 0,2 .0,4 0,6 




0,0 0,25 ^ 0,75 



FIG. 2: Numerical results for DSGE containing four and six 
sites with a "kink" centered in the middle between the sites. 
Left panels: full spectrum of a) 4 pendulums and c) 6 pendu- 
lums in the interval of kink stability; the numbers correspond 
to the mode number. Right panels show the dependencies of 
oscillations amplitude ratios of different neighbor sites corre- 
sponding to two lowest modes: b) for four pendulums V2/V1; 
d) for six pendulums Vi+i/vi with i = 1 (solid line) and i — 2 
(dashed line). The digits correspond to the mode number. 

For the symmetric modes (v n = V— n ) the first equation 
of system JHJ is to be replaced with 

(w 2 - A - cos u?)«i + \v 2 = . (9) 

For small A we use the approximate expressions for 
and in the leading approximation substitute cosmJ — > 
1 — 27r 2 A 2 , cos — * 1 for n > 1. This means that we 
effectively replaced the kink with two isotopic impurities 
located at sites 1 and -1, the "strength" of which alters 
with the change of A. Seeking the solution in the form 
v n = Ac Kn + Be~ Kn , with constant A and B, one ob- 
tains the spectral dependence for uj(k) in form J5J), where 
k = ik. Then using the consistency condition after some 
straightforward algebra we arrive at the relation which 
defines the allowed values of k for antisymmetric modes: 

(1 - 2tt 2 A) sinh[K(/V - 1)] 
+2tt 2 A sinh[/tJV] - sinh[ K (l + TV)] = . (10) 

Expanding this relation in the vicinity of k — (kN <C 
1) up to k 5 we determine the allowed values of k for 
the frequency of first antisymmetric mode, u>i, and for 
the next, third mode, W3, from the biquadratic equation: 
an 4 + bn 2 + c = 0. The expressions for the coefficients 
are: 

1(A - A d ) + 5{2N 2 + N 4 )(X - A d ) - 5A(2/V 3 + N) 

a= 60 ' 

b={\-\ d )/3 + N 2 (\-\ d )-N\, 
c = 2(A-A d ), 
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and the detachment point is Ad = 1/tt 2 . The notable fact, 
which can be extracted from Eci. (|10fl . is the following: if 
A is close to zero, A <C N" 1 , we have k 2 = — k 2 ~ iV~ 2 , 
and for this values of coupling parameter one finds the 
expressions for frequency dependencies as: 



2A(3 
2A(3 



V3)N- 2 . 
V3)N-' 2 . 



(11a) 
(lib) 



So, initially, for small A, the lowest antisymmetric mode, 
Ea. (|lla[) . goes up being inside the spectrum. However 
quite a different situation occurs if one moves inside the 
region where the inequality |A — X d \ <C iV _1 holds, i.e. 
in the close vicinity of A^. In this region one obtains 
k 2 ~ N^ 1 (here the inequality kN <C 1 is true because of 
the additional smallness provided by the factor (A — A<j)). 
Then for the frequency dependencies we have: 

u\ ps 1 - 2(A - X^N' 1 , (12a) 

lu 2 ps 1 - 4(A - X^N- 1 + 6X d N- 2 . (12b) 

The second derivative of uj\ at the point A = Xd in- 
volves the term independent on N: d 2 wi/dA 2 = — A7 1 . 
Therefore we can conclude that in the limit N — ► oo 
the splitting of this dependence from the lower bound- 
ary of the spectrum must have a parabolic form. From 
the expressions <|11I12[) it becomes evident what brings 
about the peculiarity as A approaches Xd, i.e. as the first 
antisymmetric mode detaches from the spectrum. For 
A close to zero we have dwi^/dA ~ iV~ 2 , whereas in 
the vicinity of Xd the different dependence takes place: 
dui^/dX ~ N" 1 . (Note that the sign of the first deriva- 
tive also changes.) Because of this the initial weak (~ 
N~ 2 ) growth for small values of A changes to more rapid 
(~ TV -1 ) decrease in the close vicinity of A^. The depen- 
dencies for the modulus of first derivatives at A = Ad on 
the number of sites shown in FigEl For N' = 2N > 30 
they are in a good agreement with analytical results (|12fl . 
With the increase of ./V this agreement becomes better 
inasmuch as we omitted the higher with respect to iV _1 
terms in Eas. l|12|l . In the region A < X d there must be an 
extremum point A GX t, da;i/dA|A cxt = 0, for both depen- 
dencies uji^(X), at which the monotonic growth changes 
to decreasing. These extremum points, A cx t, tend to Xd 
as N gets bigger: (Xd — A cxt ) ~ iV _1 . The dependencies 
of the value of difference, (Xd — A ox t), on the number of 
sites N are presented in the inset panel of Fig [31 For 
A > Ad the first antisymmetric mode gets into the spec- 
trum gap and becomes an internal mode. However the 
dependence for the next, third antisymmetric mode, in 
spite of its tendency to drop down, cannot cross the de- 
pendence of the preceding second (symmetric) mode and 
therefore these two dependencies get closer and closer to 
each other. 

Now consider an infinite system. We seek the solution 



in the form of localized wave, v r , 



k > 0. Where- 



upon the only condition defining the allowed values of k 
becomes as follows: 

"?=0. (13) 
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FIG. 3: The dependencies for the modulus of derivative at the 
point A = Xd on the number of sites for the first (antisymmet- 
ric) mode and for the next, third mode. The digits correspond 
to the mode number. Dashed curves for each mode show the 
analytical prediction obtained by virtue of Eqs. 1121 . The inset 
shows the dependence of deviation \d — A cx t of the extremal 
point from the detachment point as a function of number of 
sites N' = 2N. 



Expanding this relation we arrive at the dependence: 
Ki = 2A7 1 (A — Ad), and then one finds the expression 
for the frequency of the antisymmetric localized mode 
as: 



ojf = 1 . 



for A < Aw 



1-4A7 1 (A-A d ) 2 , for A > X d ■ 



(14) 



in consistency with the result for finite system. The char- 
acteristic spatial scale for the eigenfunction of this inter- 
nal mode is /i = k^ 1 . Obviously, the smaller l±, the 
better the eigenfunction is localized. We see that at the 
outset, as this mode detaches, h ~ (A — A^) -1 . 

Proceed to studying the symmetric modes. The rela- 
tion which defines the allowed values of n for those reads: 



-2(1 



(1 + 2tt 2 A) sinh[/t(AT - 1)] 
7T 2 A )sinh[K7V] + s inh[K(l + N)\ = 0. (15) 



Then we again expand this relation for kN <C 1 up to 
the fifth power and determine k from the biquadratic 
equation: dn A + gn 2 + f = 0, where 



d = 



5N(X + X d ) + WN 3 (X + X d ) - X(5N 4 + ION 2 + 1) 
60 



f = N(X + X d )~(N 2 + l/3)X, 
g=-2X. 

In the region A <C iV _1 one obtains: 



1 - 2AT 1 AiV- 1 . 



ujI ps 1 + &XN~ 



(16) 



27r 2 A-3-4Asinli 



The dependencies for symmetric modes have neither an 
inflection point nor any peculiarity at A = A^. 
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For the infinite system the relation for k is written as 
follows: 

e- K + 27r 2 A-l-4Asinh 2 ^ =0. (17) 

Then we arrive at the dependence for PN-mode for weak 
coupling in the form: 

lo 2 « 1 - 4A d - 2 A 3 , (18) 

which is in agreement with this dependence given in 
Ref.0. The localization distance now is: l — k q ~ 
^-3/2 xhus we can infer that the localization for PN- 
mode with the growth of coupling value A develops faster 
than for antisymmetric internal mode, as it was pointed 
above while considering small size systems (cf. FigEl 
right panels). 

The results for the frequency dependencies I|14I18|) 
could also be obtained with the use of Green func- 
tion (Lifshitz) technique for a linear chain with impu- 
rities, (see the monograph |9j), which was employed in 
Refs.jia. 

However, for the sake of getting simple an- 
alytical expressions one should replace the kink solution 
with the finite number of effective impurities. In the case 
of two impurities model, utilized in the current paper, one 
would obtain exactly the same results. In fact, the usage 
of the two impurities model means that we retain only 
leading terms in powers of parameter A. 

To sum up, we investigated the features of spectrum 
and kink internal modes of DSGE focusing on the effect 



of detachment of an internal mode from the quasicontin- 
uous spectrum. We note that the two-impurities model, 
used for our analysis, gives a somewhat different value 
for the detachment point (Xd = 7r~ 2 , while the numer- 
ically estimated value for DSGE is Xd ~ 0.26). This is 
explained by the fact that the value of Xd is nonzero, 
whereas we had to use the asymptotical expansions valid 
for A — > 0. This means the higher order terms with 
respect to A should always result in the corrections to 
the analytically obtained value of A^- However the two 
impurities model allowed us to obtain analytical results 
describing the spectrum peculiarities, which are in good 
agreement (qualitative and even quantitative for large N) 
with the numerical ones. For the infinite system the de- 
pendence of first antisymmetric mode detaches smoothly 
(with zero first derivative) from the band edge, but for 
the finite system it does cross the bandgap threshold. 
For extremely small values of A the first mode belongs to 
the spectrum. The initial weak growth (~ TV -2 ) of low- 
est antisymmetric modes frequencies with the value of A 
then changes to more rapid (~ iV -1 ) decreasing. The 
extremum point A cxt approaches the detachment point 
Xd as the number of sites gets bigger. At A = Xd the 
first antisymmetric mode drops into the spectrum gap. 
The higher antisymmetric mode come closer and closer 
to the preceding symmetric modes with the increase of 
coupling but this tendency weakens with the growth of 
mode number. The symmetric modes do not display any 
peculiarity at A^. 
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